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ABSTRACT 

Communications  Research  Laboratory  (CRL)  has  developed  a  GPS  codeless 
receiver  called  ”GTR-2”  for  measuring  total  electron  content  (TEC)  along 
the  line  of  sight  to  the  GPS  satellite  by  using  the  cross  correlation 
amplitude  of  the  received  P-code  signals  carried  by  L 1 ( 1 575 . 42MHz)  and  L2 
(1227.6MHz).  This  equipment  has  the  performance  of  uncertainty  in  the 
measurement  of  TEC  of  about  2  x  1016  electrons/m2  when  a  10  dBi  gain 
antenna  was  used. 

To  increase  the  measurement  performance,  CRL  is  planning  an  upper  version 
of  GTR-2  called  GTR-3  which  uses  the  phase  information  of  the  continuous 
signals  obtained  by  making  a  cross  correlation  or  multiplication  of  the 
received  LI  and  L2  P-code  signals. 

By  using  the  difference  of  these  measured  phase  values,  we  can  estimate 
the  ionospheric  delay  with  the  ambiguities  of  the  periods  of  L1+L2  and  LI 
-L2  signals.  As  the  periods  of  these  signals  are  about  3  ns  and  0.3  ns 
respectively,  then  this  method  has  the  possibility  of  TEC  measurement 
with  the  uncertainty  of  1  or  2  times  of  1015  electrons/m2  of  TEC. 

Additionally  to  the  precise  measurement  of  TEC,  this  method  has  the 
ability  of  the  precise  measurement  of  the  pseudo-range  between  the  GPS 
satellite  and  the  receiver  for  the  precise  positioning. 

1  .  INTRODUCTION 

Several  kinds  of  GPS  codeless  receiving  systems  for  the  precise  relative  positioning 
[l”4]  are  developed  and  they  use  the  dual  frequency  signals  transmitted  from  GPS  for 
the  calibration  or  correction  of  the  ionospheric  propagation  delay.  To  make  this 
correction  of  ionospheric  delay,  they  must  make  the  measurements  of  the  pseudo-range 
for  LI  and  L2  signals  individually. 

CRL  has  developed  a  TEC  measurement  system  using  the  amplitude  of  cross  correlation 
of  the  received  LI  and  L2  P-code  signalsC5].  This  method  is  very  simple,  but  its 
precision  is  principally  based  on  the  duration  of  P-code  signal,  so  it  is  difficult 
to  obtain  sub-ns  level  of  ionospheric  delay  correction.  To  realize  more  precise 
ionospheric  delay  correction,  CRL  is  planning  to  build  a  ionospheric  TEC  measurement 
equipment  by  using  the  method  described  in  this  paper. 

2.  IONOSPHERIC  GROUP  DELAY  AND  PHASE  DELAY 

In  the  meaning  of  first  order  approximation,  the  ionospheric  group  delay  At<i„,,(fc) 
and  phase  delay  AtPion(fc)  caused  on  the  electro-magnetic  waves  which  transverse  the 
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ionosphere  are  expressed  by  the  following  equations  (1)  and  (2)  respectively. 

Atglon(fJ  =  134'Ne/fq2  (ns),  (1) 

AtP,.n(fe)  =-134*N,/fe2  (ns).  (2) 

Where  Ne  is  ionospheric  total  electron  content  for  the  line  of  site  to  the  satellite 
in  electrons/m2,  and  fc  is  the  carrier  frequency  of  signal  in  Hz.  Namely  the  group 
delay  and  the  phase  delay  hav  same  magnitude,  but  the  group  delay  has  a  positive 
value  while  the  phase  delay  has  a  negative  value. 

Ne  for  the  vertical  path  changes  from  about  1  x  1018  to  1  x  1016  [electrons/m2]  at 
the  day  time  of  solar  maximum  and  at  the  night  time  of  solar  minimum,  and  it  changes 
comp  1  i cated  1  y  from  time  to  time.  Therefore  we  can  say  that  the  ionospheric  delay  is 
one  of  the  largest  error  sources  for  the  precise  positioning  and  the  time  transfer 
by  using  space  techniques. 

3  .  IONOSPHERIC  TOTAL  ELECTRON  CONTENT  MEASUREMENT  USING  GPS  SATELLITE 
3.1  Convent inal  measurement  methods 

As  the  ionospheric  delay  depends  on  the  carrier  frequency  mentioned  in  the  previous 
section,  so  the  conventional  calibration  or  correction  methods  used  in  the  precise 
relative  positioning  by  GPS  satellite  are  using  the  relative  delay  time  between  the 
LI  and  L2  signals.  To  get  the  relative  delay  time,  especially  in  the  case  of  the 
GPS  codeless  receivers,  one  makes  the  reconstruction  of  the  clock  signal  of  P-code 
and/or  carrier  signal  of  LI  and  L2  signal  individually.  Namely  one  must  have  the 
pseudo-range  measurement  circuits  for  LI  and  L2  signals  individually.  In  addition, 
in  the  case  of  carrier  reconstruction  method,  the  measurement  precision  is  very  high 
but  the  ambiguity  resolution  is  difficult  because  the  wave  length  of  the  carrier 
signal  is  very  short. 

One  of  the  other  methods  by  using  the  GPS  signals  for  the  TEC  measurement  is  the  GTR 
-2  method  which  is  developed  at  CRL  and  B1PM.  It  uses  the  cross-correlation  or 
multiplication  amplitude  of  LI  and  L2  P-code  signals.  As  the  measurement  precision 
of  GTR-2  method  is  proposed  to  duration  of  the  P-code  clock  signal  and  inversely 
proposed  to  the  square  root  of  signal-to-noise  ratio  (S/N),  one  must  use  a  high  gain 
antenna  to  improve  the  precision. 

3.2  Principle  of  GTR-3 

As  the  LI  and  L2  P-code  signals  transmitted  from  the  GPS  satellite  have  a  same  code 
pattern  for  each  GPS  satellite  and  all  of  the  carrier  frequencies  and  the  code  clock 
frequencies  are  synthesized  from  same  onboard  frequency  reference,  then  we  can 
express  the  LI  and  L2  P~code  signals  by  the  following  equations  when  they  are  trans¬ 
mitted  from  the  satellite; 

Xpi(t)  =  P(t)cos(2ffi.it),  (3) 

xP2(t)  -  P(t)cos(2rfL2t),  (3’) 

where  P(t)  is  a  binary  pseudo-random-noise  sequence  which  has  an  amplitude  of  t\, 
and  fLi  and  fL2  are  the  carrier  frequencies  of  LI  and  L2  signals. 

The  received  P-code  signals  referred  to  receiver  clock  are  denoted  by  the 
following  equations  (4)  and  (4’),  and  figure  1  shows  the  schematic  description  of 
them. 
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(9) 

(9’) 


U  =  ♦l/C2fFi) 

=  Atc- 134N*/(f Lif L2)-NiT  o-(l/4+ni)Tx, 

1 2  =  $2/(2rF2) 

=  Atc+134Ne/(fLifL2)~N2To~(l/ 4+n2)T2 , 

where 

T i  =  1/Fi  '  0.357ns, 

T2  =  1/Fa  ^  2.875ns, 

Fi  =  fn  +  fLj  =  10.23MHz  x  274  =  F0  x  137, 

F2  =  f l i ” f L2  =  10.23MHz  x  34  =  F0  x  17, 

F0  =  20.46MHz  imaximum  common  frequency  for  Fi  and  F2, 

T0  =  1  /F o  ^  48.88ns, 

Ni,N2,nl,n2  Ipositive  integer  numbers  (ambiguities), 

0  i  ni  <  137,  and  0  ^  n2  <  17. 

Fron  (9), O’),  we  can  estimate  TEC  and  psudo-range, 

Ne  *=  {(t2-ti)+(Na-Ni)To+(l/4+na)T2-(l/4+ni)Ti}fLifL3/268  (10) 

Atc  -  {t1+t2+(Ni+N2)T0+(l/4+ni)T1+(l/4+n2)T2}/2,  (11) 

with  the  ambiguities  of  nl,  n2,  Nl,  N2. 

Generally,  the  ambiguity  resolution  is  difficult  because  durations  of  L1+L2  signal 
and  L1-L2  signal  are  very  short,  but  in  the  case  of  GTR-3,  we  can  use  the  rough 
estimation  of  TEC  obtained  by  the  GTR-2  method  and  also  use  the  frequency  relation 
between  FI  and  F2,  the  maximum  common  frequency  of  them  is  20.46  MHz,  so  we  can 
solve  the  ambiguities  easily. 

3.3  ESTIMATION  OF  MEASUREMENT  PRECISION  AND  EXAMPLE  OF  RECEIVER  DESIGN 

Table  1  shows  the  estimated  precisions  of  TEC  and  pseudo-range  measurement  using 
GTR-3  method.  The  received  powers  of  the  LI  and  L2  P-code  signals  are  assumed  the 
minimum  values  of  the  GPS  specifications.  And  an  omni-directional  antenna  which 
has  a  gain  of  +3  dBi  is  assumed  for  the  receiving  antenna.  In  this  case,  the 
estimated  measurement  precisions  are  expected  to  be  about  1.2  x  1015  for  TEC  and  2.5 
cm  for  pseudo-range  measurement  respectively.  These  values  of  precision  are  not 
including  the  stabilities  of  the  satellite  reference  clock  and  the  receiver  clock. 
Figure  4  shows  the  block  diagram  based  of  the  GTR-3  method.  It  has  a  very  simple 
construction  and  similar  to  the  GTR-2,  but  to  keep  the  phase  information  of  the 
received  signals,  all  of  the  local  signals  for  frequency  conversion  in  the  receiver 
are  synchronized  to  the  receiver  reference  clock.  The  received  Li  and  L2  signals 
are  frequency  converted  to  the  IF  signals  and  made  a  multiplication  by  using  an 
analogue  multiplier.  The  output  signal  of  the  multiplier  have  the  continuous 
signals  corresponding  to  the  L1+L2  signal  and  the  L1-L2  signal  described  in  previous 
section.  As  they  have  about  +  9kHz  for  the  L1+L2  signal  and  ±  1kHz  for  the  L1-L2 
signal  of  doppler  frequency  shift  due  to  the  satellite  motion,  then  they  are  finally 
frequency  converted  to  the  signals  which  have  the  center  frequencies  of  10  kHz  to  be 
read  by  the  analogue-to-digital  converters  which  have  the  clock  rate  of  40kHz.  The 
obtained  digital  signals  are  signal-processed  by  micro-computer  to  calculate  the 
phases  or  the  zero-crossing  timing. 

The  receiving  system  based  on  this  concept  is  now  under  development  at  CRL,  and 
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yPi(t)  =  P(t-At,i)cos{2<fLi(t-AtPi)}, 
yP2(t)  =  P(t-AtB2)cos{2*fL2(t-AtP2)}, 


(4) 

(4’) 


where, 

Atg(fLn)  =  At.+p/c+AUrop+AtgionCfLn)  : pseudo-de  1  ay  for  modulation  term, 
Atp(ft,n)  =  At.+p/c+Attrop+AtpionCfLn)  ■ pseudo-de 1  ay  for  carrier  term, 

At,  : t i me  diference  between  satellite  clock  and  receiver  reference  clock, 
AttroP  -'tropospheric  propagation  delay, 
for  n^l  and  2. 

By  making  the  multiplication  of  yPi  and  yP2  using  a  multiplier  which  is  illustrated 
in  figure  2,  we  get  continuous  signals  or  de-spread  signals  at  the  output  port  of 
the  multiplier.  Equation  (5)  denotes  the  formula  expression  of  this  multiplication 
of  yPi  and  yP2.  and  figure  S  shows  the  spectrum  diagram  of  the  received  LI  and  L2  P- 
code  signals  and  the  result  of  their  multiplication. 

z(t)  =  yPi(t)yp2(t) 

=  P(t-Atgi)P(t-AtB2)cos{2*fLi(t-AtPl)}cos{2rfL2(t-Atp2)} 

=  A(t)[{cos(2*Fit-$i)+cos(2lF2t-t2)}]  (5) 


where 

A(t)  =  PCt-At.j)P(t-Atga), 

Fl  =  fLl+fL2. 

F2  =  f  L.  1  "" f  L2  , 

♦l  =  2*f LlAtpi+2tfL2Atp2.  $2  =  2f f LlAtpl~2f f L2Atp2, 

Ate  —  At.+p/c+AttroP- 

The  A(t)  in  the  equation  (5)  is  expressed  by  equation  (6)[6]. 

A(t)  =  B(Atgi-Atg2)+{compornents  of  P-code  clock) 

+{un  de-spread  part)  (6) 

The  first  term  represents  the  dc  coipornent  which  is  proporsed  to  the  relative  delay 
between  the  LI  and  L2  P-codes  due  to  the  ionosphere,  the  second  term  denotes  the 
P-code  clock  and  its  higher  order  compornents,  and  the  third  term  denotes  the 
compornent  which  can  not  be  de-spread. 

From  equation  (5)  and  (6),  we  can  obtain  the  continuous  signals  at  the  outputs  of 
the  band-pass  filters  in  the  figure  2  and  they  are  expressed  by 

zi(t)  =  B(Atgi-Atg2)cos(2gFit-$i),  (7) 

z  2  ( t )  =  B(Atgi-Atg2)cos(2*F2t-$2).  ( V ) 

The  phase  tesrms  in  equation  (7)  and  (7’)  are  discribed  as  the  followings! 

=  2fFiAtc-134N.2iFi/(fLifLi).  (8) 

♦2  =  2*F2Ate+134N,2|F2/(fLifL2).  (8’) 

Finally  we  get  the  epocks  of  the  zero-crossing  point  of  zi(t)  and  z2(t),  ti  and  t2, 
and  they  are  denoted  by 
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being  made  the  preliminary  receiving  tests,  and  the  results  will  be  shown  in  the 
other  occasion. 


4.  CONCLUT IONS 

In  this  paper  the  principle  of  the  method  for  precise  TEC  measurement  by  using  the 
signal  phase  informations  of  continuous  signal  of  L1+L2  and  L1-L2  obtained  by  making 
multiplication  of  received  11  and  L2  P-code  signal  from  GPS  satellite.  As  the  hard- 
wear  construction  is  very  simple  and  it  also  has  the  possibility  of  cm  level  of  the 
pseudo-range  measurement,  then  it  is  expected  to  be  used  not  only  for  the  precise 
TEC  measurement  but  also  for  the  precise  relative  positioning  receiver.  And  the 
ambiguity  resolutionin  the  measurement  of  TEC  and  pseudo-range  will  be  made  compara¬ 
tively  easily  by  using  rough  estimation  of  TEC  obtained  by  the  GTR-2  method  and 
also  the  frequency  relationbetween  L 1 +L2  and  L1-L2  frequencies. 
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Figure  1.  Scheiatic  expression  of  psudo-range. 
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Figure  3.  Scheiatic  specrtui  diagrai  of  received  signals  and  results 

of  their  lultipl ication. 
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Figure  4.  Block  diagrai  of  receiving  systei  which  is  under  development  at  CRL. 


Table  1.  Estiiatlon  of  perfonance  of  TEC  and  psudo-range  aeasureients. 


LI  P  code 

L2  P  code 

received  power 

>-136  dBm 

antenna  gain  (omi-directional) 

+3  dBi 

+3  dBi 

systei  noise  (No) 

-174  dBm/Hz 

-174  dBm/Hz 

carr ier-to-noise  ratio  (C/No) 
at  input  of  multiplier 

>+44  dBHz 

>+4i  dBHz 

carr ier-to-noise  ratio  (C/No) 
at  output  of  mltiplier 

>+12  dB-Hz 

bandwidth  of  zero-crossing  detection 
for  ti  and  t2 

2  Hz 

■easureient  precision  of  ti 

'•0.02  ns 

■easureient  precision  of  t2 

'0. 16  ns 

estiaation  precision  of  TEC 

■'•'1.2  x  101*  electrons/m2 

estimation  precision  of  psudo-range 

cm 

251/252 


